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Abstract. We present adaptive optics assisted spectroscopy of three quasars obtained with NACO at VLT. The high angular 
resolution achieved with the adaptive optics (~ 0.08"), joined to the diagnostic power of near-IR spectroscopy, allow us to 
investigate the properties of the innermost 100 pc of these quasars. In the quasar with the best adaptive optics correction, 
PG1 126-041, we spatially resolve the Pao- emission within the nuclear 100 pc. The comparison with higher excitation lines 
suggests that the narrow Pao- emission is due to nuclear star formation. The inferred intensity of the nuclear star formation 
(13 M yr -1 ) may account for most of the far-IR luminosity observed in this quasar. 
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1. Introduction 

The study of the stellar populations and star formation in the 
host galaxies of Active Galactic Nuclei (AGNs) is fundamen- 
tal to understand the connection between black-hole growth 
and galaxy formation. Several evidences for such a connec- 
tion have been obtained in local low-luminosity AGNs, i.e. 
Seyfert galaxies (e.g. Cid Fernandes et al. 120011 Heckman et 
al. [T997I Maiolino et al. [T9951 Oliva et al. fT999b . At higher, 
quasar-like luminosities the investigation of star formation and 
stellar population is much more difficult because the tracers of 
star formation are heavily diluted by the strong active nucleus. 
Nevertheless a few observational programs have been success- 
ful in disentangling the host galaxy and star formation from 
the quasar light. The separation of the two components has 
been achieved either through a spectroscopic decomposition 
or through high angular resolution imaging observations (e.g. 
Canalizo & Stockton 2001 Courbin et al. 2002 Dunlop et al. 
120031 Jahnke et al. l2003l Schade et al. l2000l By using such 
observating strategies, evidence both for active star formation 
and relatively quiescent hosts was found. However, all these 
studies probed the stellar population on scales of a few kpc, 
while the stellar activity within the central 100 pc, where any 
interplay with the quasar activity is expected to occur, remains 
poorly explored. 
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* Based on data obtained at the VLT through the ESO program 
71.B-0453(A). 



During recent years the performances of adaptive optics 
systems have improved significantly, achieving very high angu- 
lar resolutions (;£ 0.1") even on relatively faint sources. Such 
high angular resolution, coupled with spectroscopic capabili- 
ties of some of these systems, allows us to efficiently disen- 
tangle components associated with stellar activity even in the 
circumnuclear region of luminous quasars. In this letter we 
present NACO observations of three quasars at z ~0.06 selected 
to have an excess of far-Infrared emission. In one of the quasars 
we find evidence for star formation within the central 100 pc 
with a rate of ~ 13 M yr -1 , which may account for most of 
its far-IR emission. We assume a "standard" cosmology, with 
H = 70 km s _1 Mpc" 1 , Q m = 0.3 and O a = 0.7. 



2. Sample selection, observations and data 
reduction 

The targets were selected from the sample of Andreani et al. 
(1999 1 who provide mid-IR and far-IR data for a large num- 
ber of optically selected quasars. We selected quasars match- 
ing the following constraints: 1) z~0.06, so that the CONICA 
projected resolution on the source is better than 100 pc (in the 
K band) and, at the same time, both Paa and [SiVI]1.97/im 
emission lines are shifted into the K band; 2) far-IR luminos- 
ity Lfir > 10 44 erg s _1 , i.e. objects more IR luminous than 
the average of quasars at the same redshift and, therefore, sus- 
pected to harbor starburst activity; 3) radio quiet, to avoid syn- 
chrotron contamination to the far-IR radiation; 4) visible mag- 
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Name 


RA(J2000) 


DEC(J2000) 


z 


my 


Lfir 


Date 


seeing 


FWHM 


SR 


IRAS09 149-6206 


09 16 09.4 


-62 19 29 


0.057 


14.1 


7.9 


Apr 20, 2003 


1.0" 


0.09" (100 pc) 


5% 


PG1 126-041 


11 29 16.6 


-04 24 08 


0.060 


14.5 


3.2 


Jun 17, 2003 


0.4" 


0.07" (81 pc) 


26% 


PG2214+139 


22 17 12.3 


+ 14 14 21 


0.066 


14.5 


1.5 


Jun 18, 2003 


0.8" 


0.09" (114pc) 


11% 



Table 1. List of the bserved quasars, along with their properties and log of the observations. Lfir is in units of 10 44 erg s ; seeing 
is the median DIMM seeing during the observation; FWHM is the Full Width Half Maximum of the PSF on the Adaptive Optics 
corrected images, in arcsec and in parsec; SR is the Strehl Ratio measured on the images. 



nitude V<14.5 mag, so that the Adaptive Optics correction can 
achieve good performances with Strehl ratios of -0.2 by using 
the quasar itself as a reference to close the loop. Four quasars 
were selected and three of them, whose properties are listed 
in Tab. 1, were observed in service mode with NACO. Both 
images and spectra were obtained in the K band. Images were 
obtained in the IB_2.18 filter (to avoid saturation) with the S27 
camera, yielding a pixel scale of 0.027"/pix. Spectra were ob- 
tained with a 0.086" slit at PA=0°, with the Grism2 coupled 
with the SK order sorting filter, yielding a spectral resolution 
R = A/AA = 1400, and with the S54 camera yielding a pixel 
scale along the slit of 0.054"/pix. During each observing run 
we obtained a diffraction limited nuclear spike of ~ 0.08" (cor- 
responding to ~ 100 pc projected on these quasars). In the spec- 
troscopic mode the pixel size undersamples the PSF. The Strehl 
Ratio ranges from 5% to 26% (see Tab. 1). 

The spectra were reduced following the standard steps for 
NIR spectroscopy. The spectra obtained at different positions 
along the slit were subtracted from each other to remove the 
background, then flat fielded, aligned, co-added and calibrated. 
Atmospheric features were corrected by dividing for the spec- 
trum of a reference star. B3III and B5V stars were used for 
PG1126 and IRAS09149 respectively, and in these cases the 
spectra were then multiplied by a blackbody at 20000 K to re- 
establish the correct slope of the continuum. For PG2214 we 
used a GOV star, and the spectra were then multiplied by the so- 
lar spectrum to remove the stellar features and continuum slope 
(see Maiolino et al. 1996 for details). Images were instead re- 
duced using the eclipse software and the CONICA pipeline. 

3. Analysis and results 

In this letter we focus on the analysis of the spectra, since they 
provide the most interesting information. Of the three quasar 
spectra we obtain spatially resolved features only for the case 
with the best Adaptive Optics correction, i.e. PG1 126-041 
which was observed with a Strehl Ratio of 26%. The integrated 
spectrum is shown in Fig.^ In the following we will concen- 
trate on the analysis of this spectrum. 

Paa is clearly broad and certainly dominated by the Broad 
Line Region (BLR). It is more difficult to indentify the pres- 
ence of a narrow component associated either with a Narrow 
Line Region (NLR) or with star formation. A method to inves- 
tigate the presence of a narrow component, while exploiting 
the high angular resolution delivered by the AO, is to study the 
spatial variation of the Paa profile along the slit. The BLR is 
certainly unresolved at our resolution (R<lpc) and therefore 
the profile of the Paa, if totally attributed to the BLR, should 
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Fig. 1. Continuum subtracted spectra of PG1 126-04. The Bry* 
introduced by the division of the reference B star is marked. 



not change spatially along the slit, and should simply scale in 
intensity following the profile of the PSF. We have investigated 
the variation of the Paa profile along the slit by subtracting the 
continuum of all spectra at different locations along the slit and 
then rescaled the spectra to match the Paa flux on the central 
spectrum. In Fig.|2]we show this comparison for the spectrum 
0.054" north of the central spectrum for PG1 126: the solid line 
shows the (continuum subtracted) central spectrum, while the 
dashed line shows the rescaled off-nuclear spectrum. The lower 
panel shows the difference between the two spectra which, al- 
though noisy, shows the presence of a excess narrow compo- 
nent (FWHM ~ 200 km/ s) in the off-nuclear spectrum. This 
demonstrates the presence of a narrow component of Paa, as- 
sociated either with a NLR or with star forming regions. 

Once proved the existence of a narrow component, sepa- 
rated from the BLR, we can follow its spatial behavior along 
with other spectral features. We have fitted the Paa of PG1 126 
profile with two broad gaussians, which reproduce the broad 
(BLR) component, and a narrow component. The resulting fit 
for the central spectrum is shown in Fig. [3] In Tab. 2 we list 
the best fit parameters for the components. For what concerns 
the blend of [SiVI]1.96/im and Brd>, each of the two lines was 
fitted by using two gaussian components only (additional com- 
ponents are not required given the lower S/N for these fainter 
lines) 1 . [SiVI] in particular is well fitted by a broad and a nar- 
row component. The presence of both narrow and broad com- 
ponents of coronal lines (such as [SiVI]) is well known in 
AGNs. A component with a width intermediate between BLR 
and NLR (~ 1000-2000 km s~') is commonly observed and as- 

1 Note that the profile of the Paa cannot be used to fit the Br<5 be- 
cause different hydrogen lines originating from the BLR have often 
different profiles, like in this case. 
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Line 


Flux 


v-v sys 


FWHM 


Paa - Broad 1 
Paa - Broad 2 
Pace - Narrow 


953.3 ± 43.3 
457.1 ±41.7 
27.6 ± 7.0 


14 ±37 
37 ±32 
71 ± 32 


3003 ± 134 
1312 ± 142 
Unresolved 


BrS - Broad 
Brc5 - Narrow 


129.0 ± 10.4 
13.8 ± 6.3 


-75 ± 191 
-124 ± 123 


2277 ±258 
697 ±413 


[SiVI] - Broad 
[SiVI] - Narrow 


89.7 ± 6.8 
8.0+ 1.7 


-203 ± 333 
42 ± 89 


2797 ± 240 
360 ± 185 



Table 2. Best fit parameters for the components of the central 
spectrum of PG1 126-041, for Paa, Br6 and [SiVI]. Flux is in 
units of 10~ 16 erg s _I cm 2 ; the shift the various component 
from the systemic velocity (v- 



y 5 ) and their FWHM are in km 




2 2.02 
Wavelenght (jum) 




-+- 



-+- 



-+- 



Difference 



poor atm 




2.02 

Wavelenght (/im) 

Fig. 2. Comparison between the central spectrum and the spec- 
trum 0.054" to the north, around the Paa region. In the upper 
panel the solid gray line shows the (continuum subtracted) cen- 
tral spectrum, while the black dashed line shows the rescaled 
off-nuclear spectrum. The loer panel shows the difference be- 
tween the two spectra. 



cribed to an "Intermediate Emission Line Region" (Giannuzzo 



et al. 1995 ). The narrow component of the coronal lines comes 
instead from the "classical" NLR and extends for a few 100 pc 
(Maiolino et al. 2000 Thompson et al. 2001 1. The broad com- 
ponent requires also the profile to be asymmetric, and specif- 
ically with a blue wing. Such a blue asymmetry is also com- 
monly observed on high excitation lines from the nuclear re- 
gion of several Seyferts (e.g. Marconi et al. 1199 6 Oliva et al. 
[TO). 

The line profiles on the off-nuclear spectra were fitted by 
keeping frozen the parameters of the broad components (ex- 
cept for a scaling factor) while the narrow components were left 
free to vary. In Fig.|4]we show the intensity of the various lines 
as a function of the offset from the nucleus. The broad com- 
ponent of the Paa gives the profiles of the PSF along the slit 
(and it is also consistent with the profile of the continuum, not 



Fig. 3. Fit of the Paa profile for the central spectrum of PG1 12. 
The fit is overplotted on the original spectrum, together with 
the single gaussian components (see text). 



shown). The most interesting result, expected from the analysis 
discussed above, is that the intensity of the narrow component 
of Paa does not follow the PSF profile and shows an excess of 
emission to the North, tracing an emission line region extend- 
ing for 0.05"-0.1", or 60-120 pc, from the nucleus. 

We investigate the nature of the narrow Paa emitting region 
by comparing its flux with the observed [SiVI] emission. The 
narrow component of the latter traces the NLR as it is not pro- 
duced by star forming regions, while the narrow Paa may either 
come from star formation or from the NLR. The different trends 
of the narrow Paa and [SiVI] in Fig. 0] suggest that the two 
lines have different origins. We have further investigated this 
issue by comparing the ratio between narrow Paa and [SiVI] 
in the narrow line region of Sy2 galaxies observed by other au- 
thors (Syl cannot be used because the hydrogen lines include 
the broad component and generally a decomposition of narrow 
and broad components is not given). Paa is not observable in 
low redshift Sy2s, but it can be easily derived from Bry by us- 
ing the case B recombination ratio (which holds for the NLR): 
Bry/Paa=0.082. The additional problem is however that in Sy2 
galaxies Bry may well be contributed by star formation, since 
there are several evidences for nuclear and circumnuclear star 
formation, as discussed in Sect. 1 . For this reason we decided 
to focus on two specific templates: NGC 1068 and Circinus. 
These are very nearby Sy2's, extremely well studied, for which 
there is no evidence in the central region for significant con- 
tribution to Bry by active star formation (star formation occurs 
in rings at a radius of ~ 10" - 30"). The Paa/[SiVI] ratio in- 
ferred for the NLR is 2.2 in Circinus (Maiolino et al. 1998 1 and 
1.1 in NGC 1068 (Reunanen et al. 120031 . In PG1126 the ra- 
tio between narrow Paa/[SiVI] is 3.4 on the nucleus and 6.4 in 
the northern region at 60 pc from the nucleus. This comparison 
indicates that there is a strong excess of narrow Paa emission 
with respect to that expected from a NLR. Such a Paa excess 
is probably due to active star formation in the nuclear and cir- 
cumnuclear region, within the central ~100 pc. 

Assuming that all the narrow Paa is associated with star 
formation we can estimate the star formation rate in the central 
region of this quasar. Unfortunately, our data do not provide a 
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Fig. 4. Intensity of the various lines as a function of the offset 
from the nucleus, normalized at the flux in the central spectrum. 
The peak emission of the narrow Paa is not aligned with the 
nucleus, but shows an excess ~60 pc to the north. 



bi-dimensional information of the Paa distribution in the cen- 
tral region. However, we can roughly estimate the circumnu- 
clear integrated Paa emission by assuming that the emission is 
well represented by the average of the northern and southern 
sides in our spectrum. The inferred total narrow Paa flux esti- 
mated in this way is F tot (Pao;N all - c , w ) = 2.15 10~ 14 erg s^'ctrT 2 . 
From case B recombination (Paa/Ha = 0.107), we infer 
a total Ha luminosity due to star formation of L(Ha)sF = 
1.7 10 42 erg s _1 . By using the relation between Ha luminos- 
ity and star formation rate given in Kennicutt et al. (1998 1 
we derive a nuclear star formation rate of SFR ss 13 M 
yr . According to the relation between SFR and far-IR emis- 
sion obtained by Kennicutt et al. (1998 1, such a star forma- 
tion rate is expected to produce a far-IR luminosity of ~ 
3 10 44 erg s" 1 , which matches the luminosity observed from 
this quasar (Tab. 1). Therefore, the nuclear star formation de- 
tected in the nuclear few 100 pc may well account for most of 
the far-IR luminosity emitted by this quasar and could explain 
the excess of far-IR emission in this object (Sect. 2). A similar 
result was found at higher redshift by Alexander et al. (|2004i, 
who derived that star formation appears to dominate the bolo- 
metric output of AGNs hosted in bright SCUBA galaxies. 

Of course, the estimated nuclear star formation rate is de- 
rived from a monodimensional spectrum and, in particular, the 
inferred nuclear (R<250 pc) Paa luminosity has required a 
large aperture correction (a factor of ~3). As consequence, a 
confirmation of these findings is certainly required with inte- 
gral field spectroscopy. Another source of uncertainty is the 
possible star formation activity occurring in the host galaxy at 
radii larger than 250 pc; indeed the limited signal-to-noise ratio 
of our spectra in the outer regions can only provide a relatively 



loose upper limit of 0.18 M yr~' kpc 2 on the star formation 
rate per unit surface area. Deeper observations are required to 
further constrain the star formation in the host galaxy. 

In the other two quasars the analysis of the line profile is 
more complex due to imperfect subtraction of some deep at- 
mospheric absorption features and to a much lower Strehl ratio 
than in PG1126. There are some hints of a resolved narrow 
component of Paa, but need to be confirmed with higher qual- 
ity spectra and higher Strehl Ratios. 

4. Conclusions 

We have obtained K-band spectroscopic observations assisted 
by adaptive optics of three quasars at z~0.06 which are lumi- 
nous in the far-IR. The adaptive optics correction allows us the 
reach a diffraction limited angular resolution (~ 0.08", cor- 
responding to ~ 100 pc), with Strehl ratios ranging from 5% 
to 26%. In the quasar observed with the highest Strehl ratio 
(PG1 126-041) we spatially resolve a narrow component of 
Paa on scales of ~50-100 pc. By comparing the spatial dis- 
tribution of the narrow Paa with the distribution of the coronal 
line [SiVI] 1 .96/mi, we infer that most of the narrow Paa is due 
to star formation in the nuclear and circumnuclear region of the 
quasar. We also derive that most of the far-IR emission of this 
quasar is due to the nuclear star formation detected by us. 

These observations clearly demonstrate that adaptive optics 
assisted spectroscopy is a powerful tool to investigate the nu- 
clear region even in bright quasars. 
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